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Abstract. We present measurements of the emission from the centers of fifteen spiral galaxies in the 3 Pi — 3 Po [Cl] 
fine-structure transition at 492 GHz. Observed galaxy centers range from quiescent to starburst to active. The 
intensities of neutral carbon, the J=2-l transition of 13 CO and the J=4-3 transition of 12 CO are compared in 
matched beams. Most galaxy centers emit more strongly in [Cl] than in 13 CO, completely unlike the situation 
pertaining to Galactic molecular cloud regions. [Cl] intensities are lower than, but nevertheless comparable to 
J= 4-3 12 CO intensities, again rather different from Galactic sources. The ratio of [Cl] to 13 CO increases with 
the central [Cl] luminosity of a galaxy; it is lowest for quiescent and mild starburst centers, and highest for 
strong starburst centers and active nuclei. Comparison with radiative transfer model calculations shows that most 
observed galaxy centers have neutral carbon abundances close to, or exceeding, carbon monoxide abundances, 
rather independent from the assumed model gas parameters. The same models suggest that the emission from 
neutral carbon and carbon monoxide, if assumed to originate in the same volumes, arises from a warm and dense 
gas rather than a hot and tenuous, or a cold and very dense gas. The observed [Cl] intensities together with 
literature [CII] line and far-infrared continuum data likewise suggest that a significant fraction of the emission 
originates in medium-density gas (n = 10 3 — 10 4 cm -3 ), subjected to radiation fields of various strengths. 
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1. Introduction 


Carbon monoxide (CO), the most common molecule af¬ 
ter H 2 , is now routinely detected in external galaxies. 
However, when exposed to energetic radiation, CO is read¬ 
ily photodissociated turning atomic carbon into an impor¬ 
tant constituent of the interstellar medium. As the ion¬ 
ization potential of neutral carbon is quite close to the 
dissociation energy of CO, neutral carbon subsequently 
may be ionized rather easily. As a consequence, [Cl] emis¬ 
sion primarily arises from interface regions between zones 
emitting in [CII] and CO respectively (see e.g. Israel et 
al. 1996; Bolatto et al. 2000). It requires column densi¬ 
ties sufficiently high for shielding against ionizing radia¬ 
tion, but not so high that CO selfshielding allows most 
gas-phase carbon to be bound in molecules. In principle, 
observations of emission from CO, C° and C + provide 
significant information on the physical condition of the 
cloud complexes from which the emission arises (Israel, 
Tilanus & Baas 1998; Gerin & Phillips 2000; Israel & 
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Baas 2001). Even though the far-infrared continuum and 
the [CII] line are much more efficient coolants, the vari¬ 
ous CO and [Cl] lines are important coolants for relatively 
cool, dense molecular gas, contributing about equally to 
its cooling (Israel et al. 1995; Gerin & Phillips 2000). In 
galaxies, however, studies of the dense interstellar medium 
are complicated by the effectively very large linear ob¬ 
serving beams which incorporate whole ensembles of in¬ 
dividual, mutually different clouds. The clumpy nature 
of the interstellar medium allows UV radiation to pene¬ 
trate deeply into cloud complexes, so that the CO, [Cl] 
and [CII] emitting volumes appear to coincide when ob¬ 
served with large beamsizes. The physics and structure 
of such photon-dominated regions (PDR’s) has been re¬ 
viewed most recently by Hollenbach & Tielens (1999), 
whereas their consequent observational parameters have 
been modelled by e.g. Kaufman et al. (1999). 

[CII] emission has been observed towards numerous 
galaxies, both from airborne (the now defunct NASA 
Kuiper Airborne Observatory) and from spaceborne (the 
equally defunct Infrared Space Observatory) platforms. 
In contrast to these [CII] observations, observations of 
[Cl] emission can be performed on the ground, at least 
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Fig. 1. [Cl] and J=2-l 13 CO spectra observed towards sample galaxies. The vertical scale is T m b in K; the horizontal 
scale is velocity Vlsr in kms -1 . For all galaxies, the temperature range in [Cl] is four times that in 13 CO. 


Table 1 . Line observations log 


Galaxy 

Position 

RA(1950) Dec (1950) 

(h ms) (° ' ") 

Adopted 

Distance 

(Mpc) 

No. [Cl] 
Points 

[Cl] 

Date 

T 

sys 

(K) 

J=2-l 

Date 

13 CO 

T aya 

(K) 

J=4-3 

Date 

12 CO 
T ays 
(K) 

NGC 253 

00:45:05.7 

-25:33:38 

2.5 

20 

12/93 

3770 

12/93 

1695 

11/94 

9800 

NGC 278 

00:49:15.0 

+47:16:46 

12 

1 

07/96 

3650 

06/95 

480 

01/01 

1325 

NGC 660 

01:40:21.6 

+ 13:23:41 

13 

1 

07/96 

3065 

05/01 

350 

08/99 

3870 

Maffei 2 

02:38:08.5 

+59:23:24 

2.7 

1 

12/93 

4885 

01/96 

550 

07/96 

3700 

NGC 1068 

02:40:07.2 

-00:13:30 

14.4 

22 

07/96 

4000 

01/96 

455 

07/96 

3365 

IC 342 

03:41:36.6 

+67:56:25 

1.8 

27 

11/94 

4485 

02/89 

1440 

04/94 

2170 

M 82 

09:51:43.9 

+69:55:01 

3.25 

6 

12/93 

7200 

04/93 

335 

10/93 

9085 

NGC 3079 

09:58:35.4 

+55:55:11 

18.0 

7 

03/94 

6240 

06/95 

310 

03/94 

5510 

NGC 3628 

11:17:41.6 

+ 13:51:40 

6.7 

8 

11/94 

3450 

06/95 

325 

03/94 

2414 

NGC 4826 

12:54:17.4 

+21:57:06 

5.1 

2 

03/97 

3520 

12/93 

535 

12/93 

2045 

M 51 

13:27:45.3 

+47:27:25 

9.6 

4 

11/94 

6600 

06/95 

370 

04/96 

4065 

M 83 

13:34:11.3 

-29:36:39 

3.5 

14 

12/93 

4590 

06/95 

430 

12/93 

4360 

NGC 5713 

14:37:37.6 

-00:04:34 

21.0 

1 

02/99 

8000 

12/00 

515 



NGC 6946 

20:33:48.8 

+59:58:50 

5.5 

17 

07/96 

3970 

01/96 

530 

07/96 

2900 

NGC 7331 

22:34:46.6 

+34:09:21 

14.3 

1 

11/96 

1925 

12/97 

320 




in the 3 Pi — 3 Po [Cl] transition at 492 GHz. However, at¬ 
mospheric transparency is poor at such high frequencies 
and weather conditions need to be unusually favourable 
for observations of the often weak extragalactic [Cl] emis¬ 
sion to succeed, even at the excellent high-altitude site of 
telescopes as the JCMT and the CSO. Consequently, the 
number of published results is relatively limited. Beyond 
the Local Group, i.e. at distances larger than 1 Mpc, 


[Cl] has been mapped in bright galaxies such as IC 342 
(Blittgenbach et al. 1992), M 82 (Schilke et al. 1993; White 
et al. 1994; Stutzki et al. 1997) and NGC 253 (Israel et al. 
1995; Harrison et al. 1995), as well as M 83 (Petitpas & 
Wilson, 1998; Israel & Baas 2001) and NGC 6946 (Israel 
& Baas 2001). A survey of 13 galaxies, including limited 
radial mapping of NGC 891 and NGC 6946 was recently 
published by Gerin & Phillips (2000). In this paper, we 
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Table 2. [Cl], <7=2-1 13 CO, <7=4-3 12 CO intensities 


Galaxy 

Offset 

// 

T mb ([CI]) 

(10") 

(mK) 

7 ([Cl]) 

Center Position 
/( 12 C04 — 3) 

(K km s' 

7 ([Cl]) 
(22") 

- 1 ) 

/( 13 CO) 

Area-integrated 
[Cl] Luminosity 

(K kms -1 kpc 2 ) 

NGC 253 

0 , 0 

2615 

486A60 

1019±120 

290±45 

106±13 

14±2.2 

NGC 278 

0 , 0 

100 

7±3 

9±2 

(5±1) 

2.6A0.4 

(3.1±0.7) 

NGC 660 

0 , 0 

240 

50±8 

85±12 

(31±8) 

7.8±1.1 

(26±7) 

Mallei 2 

0 , 0 

190 

37±7 

405±50 

(20±7) 

22±4 

(0.9A0.3) 

NGC 1068 

0 , 0 

560 

109±12 

153±19 

49±9 

11±2 

53±8.5 

IC 342 

0 , 0 

1030 

54±6 

209±21 

27±7 

24±3 

1.1±0.3 

M 82 

0 , 0 

2130 

224±35 

591±95 

180±30 

60±9 

39±6.9 

NGC 3079 

0 , 0 

530 

111±18 

115±20 

(70±15) 

12±3 

143±31 

NGC 3628 

-17, +5 

265 

84±11 

110±15 

38±8 

10±2 

28±5.7 

NGC 4826 

0 , 0 

135 

11±2 


17±4 

7.8A1.6 



-20, +5 

270 

86±10 

73±9 

(47±10) 

15±2 

(13±3) 

M 51 

0 , 0 

565 

28±5 

24±4 

(13±4) 

8.4A1.9 

(8±2.7) 


-12,-12 

340 

24±5 






-12,-24 

755 

55±9 






-24, -24 

470 

55±9 





M 83 

0 , 0 

685 

83±14 

270±20 

55±8 

29±3 

3.6±0.5 

NGC 5713 

0 , 0 

<90 

2±0.4 



7.4A1.6 

1.7-2.6 

NGC 6946 

0 , 0 

465 

85±9 

216±20 

44±8 

22±3 

13±2.4 

NGC 7331 

0 , 0 

30 

2±0.3 



2.5±0.6 

0.8-1.7 


Note: Offset position (-20,+5) of NGC 4826 is actual nucleus position. 


present a similar [Cl] survey of 15 galaxies. We also ob¬ 
tained <7=2- 1 13 CO measurements for all galaxies, and 
<7=4-3 measurements for all but two. Taking overlap into 
account, this survey brings the total number of galaxies 
outside the Local Group, detected in [Cl], to 26. 


2. Observations 


All observations were carried out with the 15m James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea (Hawaii) 
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[], mostly between 1993 and 1996. The JCMT has a point¬ 
ing accuracy better than 2" r.m.s. Observations of the 
3 Pi — 3 Po [Cl] transition at v = 492.161 GHz were made 
at a resolution of 10.2" (HPBW); those of the J=2-l 
13 CO transition with a resolution of 22". The observing 
conditions can be judged from the total system tempera¬ 
tures (including sky) which are listed in Table 1 

All observations were obtained using the Digital 
Autocorrelator Spectrograph as a backend. Spectra were 
binned to various resolutions; we applied linear baseline 
corrections only and scaled the spectra to main-beam 
brightness temperatures, T m b = T^/r] m b. Line parameters 
were determined by gaussian fitting and by adding chan¬ 
nel intensities over the relevant range. For observations 
at 492 GHz we used ij m b = 0.43 up to January 1995 and 
?7mb = 0.50 for later observations. For observations at 220 
GHz, we used r ; m b = 0.69. A list of the observed galaxies 
and additional information is given in Table 1. 

3. Results 

Observational results are summarized in Table 2; a num¬ 
ber of representative profiles of both [Cl] and J=2-l 13 CO 
emission is shown in Fig. 1. For half of the galaxy sample, 
the distribution of [Cl] was mapped beyond the central po¬ 
sition (cf. Table 1). The [Cl] emission maps of NGC 6946 
and M 83 have already been published (Israel & Baas 
2001); the remaining maps will be discussed in forthcom¬ 
ing papers. In the meantime, we have used the informa¬ 
tion contained in the maps to convolve the full-resolution 
central [Cl] intensities (Column 4 of Table 2) to [Cl] inten¬ 
sities (Column 5) appropriate to the twice larger beamsize 
of the J=2-l 13 CO observations. We have likewise used 
the map information to determine total [Cl] luminosities 
(Column 7) integrated over the entire central source ex¬ 
tent. 

For all but two of the sample galaxies we have also 
maps of the J= 4-3 12 CO distribution. Again, we refer to 
published and forthcoming papers for a discussion of these 
observations. For those galaxies that were not mapped in 
[Cl], we have used the J= 4-3 12 CO maps to estimate 
the convolved [Cl] intensity and the total luminosity by 
assuming identical [Cl] and J= 4-3 12 CO distributions. 
This was shown to be the case for M 83 and NGC 6946 
(Petitpas & Wilson 1998; Israel & Baas 2001), but we 
have further verified the validity of this assumption for all 
galaxies that were mapped in both [Cl] and J= 4-3 12 CO. 
Values obtained in this way are given in parentheses in 
Table 2. 

In all galaxies mapped, the central neutral carbon peak 
is well contained within a radius R < 0.6 kpc, often as 
small as R ss 0.3 kpc. For only two galaxies (NGC 5713 

1 The James Clerk Maxwell Telescope is operated by The 
Joint Astronomy Centre on behalf of the Particle Physics 
and Astronomy Research Council of the United Kingdom, 
the Netherlands Organisation for Scientific Research, and the 
National Research Council of Canada. 


and NGC 7331) we have no information on extent. In these 
two cases we have listed [Cl] luminosities ranging from 
that observed in a single beam to that appropriate to the 
implied maximum source diameter of 1 kpc. 

The total [Cl] luminosities of the observed galaxies 
cover a large range. Quiescent galaxies (NGC 7331, IC 342, 
Maffei 2, NGC 278, NGC 5713) have modest luminosities 
« 1 < L[ci; < 5Kkms -1 kpc 2 . Galaxies with a star- 
burst nucleus (NGC 253, NGC 660, M 82, NGC 3628, 
NGC 6946) have luminosities 10 < L[ci] < 40Kkms _1 
kpc 2 . However, M 83 has only L[cil = S.OKkms^ 1 kpc 2 , 
although it is also a starburst galaxy. The highest lumi¬ 
nosities L[ci] > 50Kkms _1 kpc 2 are found in the active 
galaxies NGC 1068 and NGC 3079. 

Interestingly, the ratio of the 3 Pi — 3 Po [Cl] and J2-1 
13 CO line strengths exhibits a similar behaviour. The [Cl] 
line is stronger than J=2-l 13 CO in all galaxies except 
Maffei 2 and NGC 7331. The highest [Cl]/ 13 CO ratios of 
about five belong to the active galaxies NGC 1068 and 
NGC 3079. Generally, the 3 Pi — 3 Po [Cl] line is weaker 
than the J=4-3 12 CO line, but not by much. In NGC 278, 
NGC 3079, NGC 4826 and M 51, the two lines are roughly 
of equal strength. Only in Maffei 2 is the [Cl] line much 
weaker. 


4. Comparison of [Cl] and 13 CO intensities 

In various studies of Galactic objects, the similarity of 
3 Pi — 3 P 0 [Cl] and J=2-l 13 CO intensities and distribu¬ 
tion is commented upon. Early such studies by the CSO 
group were reviewed by Keene et al. (1996). The CSO 
mapping of the Galactic molecular cloud complexes M 17 
and Oph A show virtually identical line intensities for 
[Cl] and 13 CO throughout. This is also found for most 
of the Orion Bar and OMC-1. The densest regions of 
Orion, however, show increasingly strong 13 CO emission 
whereas 3 Pi — 3 Po [Cl] intensities level off, yielding ever 
lower [CI]/ 13 CO intensity ratios down to about 0.4. A 
similar range of ratios (0.3 - 1.1) was found by Jansen et 
al. (1996) towards the emission/reflection nebula IC 63. 
Keene et al. (1996) attributed such low values to the ef¬ 
fects of enhanced UV radiation in photon-dominated re¬ 
gions (PDR’s). This interpretation finds support in the 
results obtained by Plume et al. (1999) and Tatematsu 
et al. (1999) who used a reimaging device on the CSO to 
effectively obtain a larger beamsize suitable for large-area 
mapping. Their maps of clouds associated with the low- 
UV sources TMC-1, L 134N and IC 5146 have fairly uni¬ 
form ratios I[CI]/ 13 CO = 1.05 ± 0.15, as do the translu¬ 
cent regions of the dark cloud L 183 observed by Stark et 
al. (1996). In contrast, maps of the molecular clouds asso¬ 
ciated with the high-UV sources W3, NGC 2024, S140 and 
Cep A yield I[CI]/ 13 CO ratios of about 0.5 for the bulk 
of the clouds. However, even here intensity ratios of about 
unity are found once again at cloud edges. The distribu¬ 
tion of cloud-edge ratios even has a tail reaching a value 
of four. Only in a few globules associated with the Helix 
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Fig. 2. Left: [Cl]/(J = 2—1 13 CO) ratios versus area-integrated [Cl] luminosity L[CI\. Right: [CI ]/(J = 4 — 3 12 CO) 
ratios versus L[CI], The J([CI])//( 13 CO) ratio appears to be a well-defined function of log L([CI]); /([CI])//(4-3 12 CO) 
ratio is not. Galactic sources (not shown) would all be crowded together in the lower left corner. 


planetary nebula (Young et al. 1997) have such relatively 
high ratios of 3-5 also been found. 

Our own data on star formation regions corroborate 
this: towards the Galactic HII regions W 58 and ON-1 
(unpublished) as well as the LMC regions N 159 and N 160 
(Bolatto et al. 2000) we find intensity ratios [CI]/J=2-1 = 
0.2-0.6 for the PDR zones associated with these starform¬ 
ing regions. The two objects (W-58C and N 159-South) 
where star formation has not yet progressed to a domi¬ 
nating stage, in contrast, yield ratios of about unity. 

As Fig. 2 shows, only a few of the observed galaxy 
centers obey the same linear correlations between [Cl] and 
13 CO that characterize Galactic clouds. Fully two thirds 
of the galaxy sample has 3 Pi — 3 Po [CI]/J = 2 — 1 13 CO 
ratios well in excess of unity; the galaxies thus have much 
stronger [Cl] emission than the 13 CO intensity and the 
Galactic results would lead us to expect. 

The galaxy sample, observed at 15" resolution, dis¬ 
cussed by Gerin & Phillips (2000) has only a little overlap 
with ours, but it shows the same effect: more than two 
thirds of the positions plotted in their Fig. 7 has a ra¬ 
tio [CI]/ 13 CO > 2. For the galaxy NGC 891, Gerin & 
Phillips (2000) observed various positions along the ma¬ 
jor axis, in addition to the central region. At the distance 
of the galaxy, their beam corresponds to a linear size of 
0.5 kpc. Whereas the [Cl] intensity generally drops with 
increasing radius, the [Cl ]/ 13 CO intensity ratio increases, 
or more specifically, this ratio increases from about 2 at 
the central positions brightest in [Cl] to about 4-6 at the 
disk positions weakest in [Cl]. 

Qualitatively, low ratios are expected from regions 
which have low neutral carbon abundances. Low neutral 
carbon abundances will be found in high-UV environments 
where neutral carbon will become ionized, and in environ¬ 
ments with high gas densities and column densities. Here, 


neutral carbon disappears because of the concomitant 
higher CO formation rates at high densities and the much 
more efficient CO (self)shielding at high column densi¬ 
ties. Because of its lower abundance, 13 CO requires larger 
column densities for for efficient shielding. Conversely, in 
environments characterized by low gas column densities 
and mild UV radiation fields, such as found in translu¬ 
cent clouds and at cloud edges, CO will be mostly dissoci¬ 
ated, and most gas-phase carbon may be neutral atomic. 
The resultant relatively high neutral carbon abundance 
will then explain high [Cl] / 13 CO intensity ratios. In this 
framework, our observations and those obtained by Gerin 
& Phillips (2000) imply that most of the emission from 
galaxy centers does not come from very dense, starform¬ 
ing molecular cloud cores. 


5. [Cl]/CO modelling 

In a number of studies (e.g. Schilke et al. 1993, Tauber 
et al. 1995, Petitpas & Wilson 1998) column densities 
have been calculated assuming [Cl] and CO emission to 
occur under optically thin LTE conditions in the high- 
temperature limit. From eqns. (1) and (2) by Tauber et 
al. (1995), it follows that: 

I([CI))/I( 13 CO) = 0.0063 A 32 /(Tex) N[CI])/N(CO) 


where 

f(T ex ) = T ex /(e 7 ^ + 3e -16 ' 6 / Tex + 5 e - 55 - 5 / T -) 

while we assume an isotopic abundance A\ 3 = 

[ 12 CO]/[ 13 CO] = 40 (cf. Mauersberger & Henkel 1993; 
Henkel et al. 1998). In Fig. 3 we have marked by a 
dashed line the expected line intensity ratios for the 
case N[CI])/N(CO) = 1. In order to verify the correct¬ 
ness of the assumptions of low optical depth and LTE, 
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Fig. 3. Model line ratio I(CI)/I(J = 2—1 13 CO) versus 
Ikin resulting for equal column densities N(CI)/dV = 
N(CO)/dV (solid lines). Curves for a range of total gas 
densities are marked in particles per cc. Dashed line: the 
ratio I(CI)/I(J = 2 — 1 13 CO) versus T ex , likewise requir¬ 
ing equal column densities N(CI)/dV = N(CO)/dV but 
assuming optically thin LTE conditions. 

we have used the Leiden radiative transfer models de¬ 
scribed by Jansen (1995) and Jansen et al. (1994) to cal¬ 
culate for gas volume densities ranging from 10 2 cm' 3 to 
10 6 cm' 3 the [CI]/ 13 CO line intensity ratio correspond¬ 
ing to unit column density ratios. The calculation was 
performed for a representative value of the column den¬ 
sity, N[CI]/dV = N{CO)/dV = 1 x 10 17 cm” 2 kins" 1 ” 1 
(cf. Israel & Baas 2001). Note that in this calculation, 
the temperature parameter is the kinetic temperature T^ n 
instead of the excitation temperature T ex . For neutral car¬ 
bon, the two are not very different under the conditions 
considered, but for 13 CO the excitation temperature is 
generally much lower than the kinetic temperature over 
most of the relevant range. Under LTE conditions, the 
[Cl]/ 13 CO ratio continuously increases with temperature 
T ex . In contrast, the radiative transfer calculation shows 
that this ratio is only weakly dependent on temperature 
above Tki n « 30 K, and in fact decreases slowly with in¬ 
creasing temperature for densities up to n/EU) « 10 4 . As 
Fig. 3 illustrates, the assumption of comparable excitation 
temperatures for 13 CO and [Cl] is valid only for very high 
densities n/LU) > 10 6 cm' 3 which are unlikely to apply 
to our observed sample. 

6. [Cl] and CO column densities 

To further investigate the physical conditions character¬ 
izing the central gas clouds that give rise to the ob¬ 
served emission, we have plotted for our galaxy sam¬ 
ple the [CI]/J = 2 — 1 13 CO line intensity ratio ver¬ 


sus the [CI]/J = 4 — 3 12 CO line ratio. For comparison 
purposes, we have added points corresponding to a few 
Galactic starforming regions (White & Sandell 1995; Israel 
& Baas, unpublished), the N159/N160 starforming com¬ 
plex in the Large Magellanic Cloud (Bolatto et ah, 2000), 
and the Milky Way Center (Fixsen et al. 1999). As the 
latter do not list 13 CO intensities, we have assumed a 
J = 2 — 1 12 CO/ 13 CO intensity ratio of 8.5, which is the 
mean value we find for the galaxies observed by us (Israel 
& Baas 1999, 2001, as well as papers in preparation). 

To put the observed points in context, we have used 
the Leiden radiative transfer models to calculate the same 
line intensity ratio in a grid with gas densities in the range 
n = 500 10 000 cm' 3 , kinetic temperatures in the range 

Tkin = 10 150 K and CO column densities 7V(CO)/dV 

= 0.3, 1.0 and 3.0 xlO 17 cm -2 /kms -1 respectively. We 
considered ./V([CI])/./V(CO) abundance ratios of 0.1, 0.3, 
1.0 and 3.0 respectively. The results are shown in Fig. 4, 
always assuming an isotopic ratio [ 12 CO]/[ 13 CO] = 40. 
Small variations in the assumed isotopic ratio lead to small 
shifts in the various curves depicted in Fig. 4, mostly along 
lines of constant temperature. 

It is immediately clear from Fig. 4 that the predicted 
[CI]/ 13 CO intensity ratio is roughly proportional to the 
A r ([CI])/A r (CO) abundance ratio at any given gas-density. 
Variation of the actual CO column density by over an 
order of magnitude or variation of the gas kinetic tem¬ 
perature has very little effect on the line intensity ra¬ 
tio except at the highest densities and column densities 
where saturation effects caused by high optical depths 
become dominant. At given column densities, however, 
the [CI]/ 13 CO intensity ratio does depend on the gas- 
density and is roughly inversely proportional to \Jn. The 
[CI]/J = 4 — 3 12 CO intensity ratio strongly varies as a 
function of gas kinetic temperature and density, as well as 
column density. 

Further inspection of Fig. 4 shows that the starform¬ 
ing regions in the Milky Way and the LMC are found 
distributed along curves that mark neutral carbon ver¬ 
sus CO abundances N(C°)/NCO) « 0.1-0.3. The galaxy 
center ratios, in contrast, mostly seem to imply signifi¬ 
cantly higher neutral carbon abundances. Only the point 
representing the quiescent bulge of NGC 7331 appears 
to be associated with an equally low carbon abundance. 
Depending on the assumed value of the total gas den¬ 
sity, centers of quiescent galaxies are associated with car¬ 
bon abundances N(C°)/NCO) ~ 0.3 (n = 500 cm -3 ) to 
1.0 (?i = 10 4 cm -3 ). This is consistent with earlier deter¬ 
minations such as N(C°)/NCO) « 0.8 (-0.4, +0.7) for 
the Milky Way (Serabyn et al. 1994). In contrast, active 
galaxies have C° column densities well exceeding CO col¬ 
umn densities independent of the gas parameters assumed. 
The diagonal distribution of galaxy points roughly follows 
lines of constant kinetic temperature. The corresponding 
temperature value varies as a function of density n and 
column density (N): Tki n > 150 K for n = 500cm -3 , 
whereas X^m = 30 — 60 K for n = 0.3 — 1.0 x 10 4 cm -3 , 
N < 10 17 cm -2 /kms -1 . Only the high-density models im- 
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Fig. 4. Observed Line intensity ratios [Cl]/ 13 CO versus [CI]/CO 4-3 compared to radiative-transfer model ratios at 
selected gas densities, for a range of temperatures and column densities. All model calculations assume an isotopic 
ratio [ 12 CO]/[ 13 CO] = 40. Galaxy centers are marked by filled hexagons, LMC star formation regions N 159 and 
N 160 by open hexagons, Galactic star formation regions W 58 and ON-1 by open triangles and the Milky Way Center 
by a cross. Lines indicate families of column density ratios N[CI]/N(CO) = 0.1, 0.3, 1 and 3. Within each family, 
dotted line corresponds to _/V(CO)/dF = 3 x 10 16 cm -2 /kms -1 , solid line to !V(CO)/dF = 1 x 10 17 cm -2 /kms -1 
and dashed lines to -/V(CO)/dF = 3 x 10 17 cm -2 /kms -1 . Temperatures of (from left to right 150, 100, 60, 30, 20 and 
10 K are marked by small open circles on each curve. 


ply a kinetic temperature range covering the fairly narrow 
dust temperature range 33 K < Td < 52 K character¬ 
izing these galaxy centers (Smith & Harvey 1996). This 
can be taken as a suggestion that at least the molecu¬ 
lar carbon monoxide emission from galaxy centers arises 
mostly from warm, dense gas as opposed to either hot, 
tenuous gas or cold, very dense gas. Possible exceptions 
to this are NGC 278 and in particular NGC 7331, M 51 
and NGC 4826 which occupy positions in the diagrams 
of Fig. 4 suggesting low temperatures Tki n = 10 — 20 K 
and consistent with the full density range including the 
highest densities. 

For M 82, Stutzki et al. (1997) estimated from the 
directly observed 3 P 2 — 3 Pi [Cl]/ 3 Pi — 3 Po [Cl] line ratio 
a density n > 10 4 cm -3 and a temperature T = 50 100 

K. This is in very good agreement with our estimates. 
However, the /([CI])//( 13 CO) ratio of three suggests an 
abundance A r [CI]/A r (CO) = 2, i.e. four times higher than 
estimated by Stutzki et al. (1997), although not ruled out 
by their results - see also Scliilke et al. (1993). 

7. [Cl], [CM] and FIR intensities 

Parameters indicating a warm and dense gas is more or 
less what is expected if the emission arises from photon- 


dominated region (PDR’s - see e.g. Kaufman et al. 2000). 
In order to explore this possibility, we have produced 
Fig. 5, which presents a comparison of the 3 Pi — 3 Po [Cl] 
and 2 P 3/2 — 2 P 1/2 [CII] line and far-infrared continuum 
intensities; this Fig. 5 is directly comparable to Fig. 8 by 
Gerin & Phillips (2000) who performed a similar compar¬ 
ison. The majority of the [CII] intensities were taken from 
Stacey et al. (1991), as were the [CII]/FIR ratios. The 
Milky Way data were taken from Fixsen et al. (1999), 
those for the LMC from Israel et al. (1996) and Bolatto 
et al. (2000), and those for NGC 278 from Kaufman et al. 
(2000), who use this galaxy as an example for the appli¬ 
cation of their PDR model calculations. The [CII] mea¬ 
surements by Stacey et al. (1991) were obtained in a 55" 
aperture, i.e. much larger than the [Cl] beam. We have 
therefore used the area-integrated [Cl] intensity, and as¬ 
sumed that all [CII] emission from the central source was 
contained within the 55" aperture. As the extent of the 
central [Cl] source is always less than this, that seems 
to be a reasonable assumption. If the [CII] emission is 
more extended than that of [Cl], the relevant [CII]/[Cl] ra¬ 
tio in Fig. 5 should be lowered correspondingly. However, 
we do not believe such a correction will change the pic¬ 
ture significantly. The [Cl]/FIR ratio was obtained from 
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Fig. 5. Ionized carbon [CII] line to far-infrared continuum 
(FIR) ratio (top) and [CII]/[Cl] line ratio (bottom) as a 
function of neutral carbon [Cl] line to FIR ratio. In both 
diagrams, the position of the Milky Way center is marked 
by a cross and the positions of Magellanic Cloud objects 
by open hexagons. 


the [CII]/[Cl] and [CII]/FIR ratio; thus [Cl]/FIR may be 
somewhat higher than plotted. 

In Fig. 5, there is no longer a clear distinction be¬ 
tween various types of objects such as we found in Fig. 4. 
Rather, the [CII], [Cl] and FIR intensities define a distri¬ 
bution in which LMC star formation regions, low-activity 
galaxy centers and high-activity galaxy centers are all in¬ 
termingled. Nevertheless, the result shown in Fig. 5 bears 
a close resemblance to the the results obtained by Gerin & 
Phillips (2000). As the [Cl]/FIR ratio increases, so does 
[CII]/FIR, but not the [CII]/[Cl] ratio which decreases 
with increasing [Cl]/FIR. Qualitatively, this may be ex¬ 
plained by PDR process along the line discussed by Gerin 
& Phillips (2000). The horizontal location of the points 
in the two diagrams suggest fairly intense PDR radiation 
fields of about 300 to 1000 times the average UV radi¬ 
ation held in the Solar Neighbourhood. For the merger 
galaxy NGC 660 we have only upper limits (log [CII]/FIR 
< -3.2, log [CII]/[Cl] < 2.1) which place this galaxy in 
the same diagram positions as the ultraluminous merg¬ 
ers Arp 220 and Mrk 231 observed by Gerin & Phillips, 
which correspond to strong radiation Helds and very high 
gas densities. 

The PDR models shown in Fig. 8 by Gerin & Phillips 
provide the highest [CII]/FIR ratios for model gas densi¬ 
ties n = 10 3 — 10 4 cm -3 . Fully half of our observed ratios 
are well above the corresponding curves, although they are 
not quite as high as the ratios observed for the three LMC 


starforming regions. Note that (the limits to) the quies¬ 
cent cloud LMC N159-S in Fig. 5 likewise suggest high 
densities but only weak radiation Helds, in good agree¬ 
ment with Israel et al. (1996) and Bolatto et al. (2000). 
For many of the galaxies and for the LMC starforming 
regions, the ratio of [CII] to [Cl] intensities appears to be 
higher than predicted by the PDR models considered. For 
the LMC objects, this was already noted and discussed by 
Israel et al. (1996). They explain this situation by an in¬ 
creased mean free pathlength of energetic UV photons due 
the lower metallicity of the LMC. However, galaxy cen¬ 
ters have, if anything, a higher metallicity (see Zaritsky et 
al. 1994). A possible explanation for the apparently sim¬ 
ilar behaviour of many galaxy centers may be a greater 
degree of filamentary or cirruslike structure. In spite of 
high metallicities, this would still allow for an effectively 
increased penetration depth of UV photons. If enhanced 
exposure results in a significantly larger fraction of carbon 
atoms becoming ionized, it would explain higher [CII] to 
[Cl] emission ratios. 

So far we have assumed homogeneous media, i.e. we 
have assumed all CO, [Cl], [CII] and FIR emission to orig¬ 
inate from the same volume. This provides in a relatively 
simple manner good estimates of the physical paramters 
characterizing the interstellar medium in the observed 
galaxy centers. 

The LMC observations, which correspond to linear res¬ 
olutions one to two orders of magnitude higher than the 
galaxy center observations, illustrate that homogeneity is 
not the case. The maps shown by Israel et al. (1996) and 
Bolatto et al. (2000) show that different locations in the 
observed regions are characterized by strongly different 
emission ratios indicating domination by different ISM 
phases (i.e. neutral atomic, ionized, molecular). A sim¬ 
ilar state of affairs applies to the Galactic Center region 
(Dahmen et al. 1998). Ideally, the observations should thus 
be modelled by physical parameters varying as a function 
of location in a complex geometry. Practically, we may 
approach reality by assuming the presence of a limited 
number of distinct gas components. The analysis of mul¬ 
titransition 12 CO, 13 CO and [Cl] observations of galaxy 
centers such as those of NGC 7331, M 83 and NGC 6946 
(Israel & Baas 1999; 2001) suggests that, within the ob¬ 
servational errors, good fits to the data can be obtained 
by modelling with only two components: one being dense 
and relatively cool, the other being relatively tenuous and 
warm. 

The galaxy points in Fig. 5 can all be reproduced by 
assuming appropriate combinations of dense and cold gas 
(having high [Cl]/FIR and [CII]/FIR ratios) with strongly 
irradiated gas of lower density (low [Cl]/FIR and high 
[CII]/[Cl] ratios). The distribution of points in Fig. 5 
would thus not directly indicate the physical condition 
of the radiating gas, but rather the relative filling fac¬ 
tors of the two components. A similar argument can be 
made to solve the apparent discrepancy between the rel¬ 
atively high kinetic temperatures suggested by Fig. 4 and 
the more modest dust temperatures referred to before. In 
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the same vein, a multi-component solution requires some¬ 
what lower beam-averaged [CI]/CO abundances than sug¬ 
gested by Fig. 4. The dataset presented in this paper is, 
however, not sufficiently detailed to warrant a more quan¬ 
titative analysis such as we have presented for NGC 7331, 
M 83 and NGC 6946 (Israel & Baas 1999, 2001), and will 
present for half a dozen more in forthcoming papers. 

8. Conclusions 

1. We have measured the emission from the 492 GHz line 
corresponding to the 3 Pi — 3 Po [Cl] transition in the 
centers of fifteen nearby spiral galaxies. In the same 
galaxies, we have also measured <7 = 4 — 3 12 CO and 
J = 2 — 1 13 CO intensities for comparison with the 
3 Pi — 3 Po [Cl] line within the framework of radiative 
transfer models. 

2. Rather unlike Galactic sources, the external galaxy 
centers have [Cl] intensities generally exceeding <7=2-1 
13 CO intensities, and in a number of cases approaching 
J= 4-3 12 CO intensities. 

3. The highest area-integrated (i.e. total central) [Cl] lu¬ 
minosities are found in the active galaxies NGC 1068 
and NGC 3079. Slightly lower luminosities occur in 
strong starburst nuclei, such as those of NGC 3628 
and NGC 6946. Quiescent and weak-starburst nuclei 
have [Cl] luminosities an order of magnitude lower. 

4. The observed [Cl], 13 CO and 12 CO line ratios, inter¬ 
preted within the context of radiative transfer models, 
suggest that the bulk of the observed emission arises 
in gas with densities n > 3000 cm -3 and kinetic tem¬ 
peratures 7ki n « 30 — 60 K. Depending on the actual 
density n, most galaxy centers should have abundances 
1V([CI])/1V(C0) = 1-3, i.e. [Cl] columns just exceed¬ 
ing those of CO. Only relatively quiescent galaxy cen¬ 
ters such as those of Maffei 2, IC 342 and NGC 7331 
have abundances 1V([CI])/1V(C0) « 0.3 — 1.0 and are 
dominated by CO just as the comparison starforming 
regions in the Milky Way and the LMC. 

5. The observed [Cl] intensities, together with literature 
[CII] line and far-infrared continuum data, likewise 
suggest that a significant fraction of the emission orig¬ 
inates in medium-density gas (n = 10 3 — 10 4 cm~ 3 ), 
subjected to radiation fields of various strengths rang¬ 
ing from a few times to several thousand times the 
local Galactic radiation field. 
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